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Forty years ago, Gardner and Ashby^, May^^^and Daniel and Mackay^ studied network sta- 
bility by analysing large networks in which species interact at random. They found a negative 
correlation between system stability and system complexity. However, their works did not ad- 
dress the stability criterion of networks with equal connectance or connectivity, nor they dis- 
cussed networks with specified topology. Reliability of network is another important quan- 
tity that can be either investigated numerically or using the theory of signature^ ^, though 
in general the comparison of the reliability of two large networks with same connectance is 
difficult. Here we investigate the stability and reliability of a class of simple networks with 
a specific topology, in the form of twisted-ring. We show how that the eigenvalue spectrum 
of these simple networks can be used to classify their stability and reliability. In agreement 
with previous results, we find that the most important factor in determining the stability and 
reliability is the largest eigenvalue, which reflects the connectance of the network. However, 
for different networks with the same connectance, their eigenvalue gaps, which reflect the 
level of symmetry of the networks, provide a guide for the ordering of stability and relia- 
bility. These conclusions on the twisted ring networks is conjectured to be general, and has 
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been verified numerically by Erdos-Renyi networki^3Hl The importance of this theoretical 
tool fies in comparing the stability and reliability of two large networks designed with the 
same connectance, when numerical methods could be time-consuming. 

Scientists studying complex systems have found the language of network very successful in 
characterizing their structures. For many evolving biological, social, and ecological systenJl^^, 
the issue of their stability and reliability are of great importance as the dynamics is affected by 
the changes in the corresponding network, while the changes in the networks also affect the in- 
teractions between the components of the complex systems. Stability of a network concerns the 
way a small disturbance at one node spreads over the entire network. The problem of stability of 
a random network has been addressed in the 1970's, by Gardner and Ashb)^^, Ma)^!], and Daniel 
and Macka)^ and more recently by ecologists AUesina and Tan^^^ and by Jirsa and Din^ for time 
delayed systems. In comparison, the reliability problem is simpler as it does not involve dynam- 
ics, though it is still of great importance in the applications to communication system^, sewage 
systems and oil and gas lines^'. One way to increase reliability with extra cost is to introduce 
redundancy at the level of the components and connections. A balance between cost and system 
reliability can be achieved by design, through the search for an optimal topology of the network. 
Since both the stability and the reliability issues are intrinsically comparative, we need a tool to 
compare the stability and reliability of two networks. Here we find a simple mathematical tool to 
distinguish two networks through their eigenvalue spectrum. 

Real complex networks often are large, thus we illustrate our mathematical results using a 
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Figure 1: Two examples of a "twisted ring" R{M, N) topology of = 6 cells and M twist, with 
12 nodes (dark disks) and 24 links (lines). The dotted lines refer to the periodic boundary condition 
to form the ring by identifying the nodes at the end of the band. Which topology is more stable 
(reliable)? 

simple physical model in the form of twisted rings i?(M, N) with N unit cells and M twists (see 
supplementary information and Fig.l). Although these rings are sufficiently simple for mathemat- 
ical analysis, yet they reveal features that can be extended to both the problem of reliability and of 
stability for general networks, including random ones. 

Most of the studies of the stability of networks begins with May's celebrated resultd^E that 
address the stability condition of a large ecosystem with N interacting species, described by an 
interaction matrix simplified by a random matrix with zero mean and variance The system is 
stable if a < etc = 1/ VNC and unstable otherwise. Here C is the connectance of the network. It 
is the probability that two given nodes are connected. Hogg, Huberman, and McGlade extended 
May's analysis to a network of interacting species participating in a non-zero sum game with 
cooperation and competitiorP^I. As these works do not address the more realistic case of a network 
with given topology, for example the rings in Fig.l, we find it necessary to extend their works 
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to one described by a well-defined adjacency matrix with element aij = if there is no link 
connecting node i and node j. Values at nonzero entries are then taken from a normal distribution. 
The networks of the complex systems we considered here are more general, while previous studies 
correspond to the special case of a fully connected network with random weights. 

We now analyse the stability and reliability of the twisted rings shown in Fig.l using spectral 
analysis of their adjacency matrix. In Fig.l, the ring is formed by connecting the bands' head and 
tail, with 2N nodes and 4N links. The symmetry group of the ring networks is given by C2 x D^, 
where C2 is the cyclic group of order 2, and Dn is the dihedral group of order 2N. Note that this 
result is independent of the number of twists M of the ring. For a twisted ring, R{M, N), M > 
the eigenvalues of their adjacency matrices are given by 

When is odd, 

,227r, ,2zM7r, 
2cos( — ) ±2cos(^^),0 < t<N-l. (1) 

When N is even. 



2 cos(^) ± 2 cos(^), 0<i<N-l,M^N/2 
2 cos(^) ± (-1)*, <i < N -1,M = N/2 



We can also find the stability landscape of the 8-unit rings numerically (Fig. 2). 



(2) 



It can be seen that the stability of the system goes through a sudden drop in the small range 
of variance value, 0.05 < a < 0.20, which is similar to results by May^. Let's introduce an 
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Figure 2: The stability landscape of 8-unit rings. The curves labelled in the first five symbols in 
the legends are the stability curves of the 8-unit rings, while the curves labelled with the last two 
symbols are the stability curves of a fully-connected network with the same size (16 nodes), but 
with interaction strength reduced to the connectance factor of 3/16 = 0.19 (*) and 4/16=0.25 (-I-) 
respectively. 

"effective interaction strength" between the various species in a ring-type network as compared 
to a random matrix of the same size, with ofe// = ctC. This effectively squeezes the Gaussian 
distribution of the random matrix with a reduced interaction strength by the factor C < 1, due to 
the entries element fljj = in the adjacency matrix of the ring. For our ring network R{M, N), 
C = 3/{2N) when M = or M = N/2, and C = 4/(2iV), when < M < N/2. Our 
numerical results on the rings are similar to the original results for the stability behaviour of a 
fully-connected network with the appropriate effective interaction strength. The critical value for 
interaction strength is simply reduced by the connectivity of the rings. We can further extend our 
analysis for rings with same connectance, namely for M = 1,2,3 when = 8. Numerically 
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we find the stability for rings R{M, N) with different twist number M, but same A^(= 8) has the 



ranking to be 



R{2, 8) < 8) ^ R{3, 8) < R{0, 8) ^ i?(4, 8) (3) 
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Figure 3: Stability of 8-unit rings for M = 2,3 normalized by the stability for i?(l,8) versus 



interaction strength a. 



Observe that the two systems with highest ranking in stability, i?(4, 8) and -R(0, 8), are sys- 



tems with the lowest index (defined as the largest eigenvalue in the eigenvalue spectrum). For a 



regular graph, such as our rings, the lower the index, the lower its connectance because the index 



of a regular graph is the degree of its nodes. Here the index of i?(4, 8) and i?(0, 8) is 3 while the 



index of R{M, 8) : M = 1, 2, 3 is 4. Therefore, we expect that the stability of the ring system 



has a negative correlation with its connectance, which agrees with the conclusion of Gardner and 



Ashby Next, we like to understand the stability ranking of -R(M, 8) for M = 1, 2, 3. These three 



rings all have index 4 and they all have the same number of nodes and links. Can we predict which 



one is more stable? We compute the spectral gap of these rings and find that the gap G{M, 8) are 
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0.5858, 1.1716, 2, 1.1716, 0.5858 for M = 0, 1,2,3,4. Here the gap is the difference between 
the largest eigenvalue and the second largest eigenvalue of the adjacency matrix. The ordering 
of the size of the gap is G{2, 8) > ^(1, 8) = ^(3, 8) > G(0, 8) = G(4, 8). We see that this is 
in complete anti correlation with the ranking of their stability obtained numerically, in agreement 
with the ranking shown in Fig. 3. 

Next let's consider reliability, which is somewhat opposite to that of stability. This is intuitive 
as a fully connected network should be most reliable, but also most unstable in the sense that a 
small disturbance in any node will quickly spread to the entire network. To address reliability, 
we follow the well-established framework in signature theor)!^. Let's assume that every link of 
the network has a certain robustness p, so that the probability of failure is 1 — p. The success 
of the whole system is defined as the path-connectedness of the network. Thus, the analysis of 
network reliability is easier than that of stability in that it does not involve the dynamics between 
the various nodes of the networks. For our rings, we can numerically compute the probability of 
system working as a function of p, which is shown in Fig.4 and 5. 

The reliability ranking of all the 8-unit ring systems is given by 

R{2, 8) ^ R{1, 8) ^ R{3, 8) > R{0, 8) ^ /?(4, 8) (4) 

The most reliable system is i?(2,8), which is the least stable, while the least reliable sys- 
tem are R{0, 8) and i?(4, 8), which are the most stable. The ranking of reliability is in complete 
correlation with the gap ordering. 
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Figure 4: Reliability landscape of 8-unit rings Figure 5: Reliability of R{M, 8) normalized by 
versus the robustness p of the link for R{M, 8). the reliabihty of 8) for M = 2, 3. 

The dependence of stability and reliability on connectance is clear, but then why do we have 
such a mysterious correlation to the ordering of the spectral gap of a network? This question can 
be answered by Cheeger's Inequality, which states that the spectral gap of a network is directly 
related to the distribution of its nodes and links and hence its symmetry: 



where G is a finite, connected, (i-regular graph, A is its second largest eigenvalue, h{G) is the 
Cheeger constant of G. Graphs with larger Cheeger constant have higher expansion and vice 
versa. This inequality was proved by Dodziuk^^, and independently by Alon-Milman^"*, and 
Alon^^. Hence, for a network with larger spectral gap (= (d — A)), it would have a higher ex- 
pansion, and tends to have a higher symmetry, in the sense that the distribution of its nodes and 



d-X 



< h{G) < v/2d(I^A) 



(5) 
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links are more homogeneous so that there are less bottlenecks in the network; whereas for a net- 
work with smaller gap would have a lower expansion and lower symmetry, corresponding to a less 
homogeneous network with bottlenecks. Therefore, the observed dependence of network stabil- 
ity and reliability on the largest eigenvalue and the spectral gap is actually a manifestation of the 
dependence on the connectance and symmetry of the network respectively. When comparing two 
networks for their stability and reliability, the first thing to check is the connectance. Networks 
with higher connectance would have a lower stability and higher reliability. If they have the same 
connectance, then we have to check for the second order effect, namely the symmetry of the net- 
works. The network with a larger gap corresponds to one with higher symmetry, and will have 
lower stability and higher reliability. These observations have been verified for the twisted rings, 
Mobius strip, as well as random networks. We conjecture that the above guideline in comparing the 
stability and reliability of two networks is applicable in general (see supplementary information). 
This theoretical tool is quite powerful, especially for comparing two large networks designed with 
same connectance, when other analytical or numerical methods are much more time-consuming. 
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Supplementary Information 



Eigenvalue Spectra of Twisted Rings We introduce some notations. We call an n x n matrix 



circulant, and denote it by drc(ai, 02, ■■ ■ , a„) if it is of the form 



circ[ai, 02, ■ ■ ■ ? 



ai 02 0,3 ' ' ' O'n 



a„ cii 0,2 



0,2 O3 Cli 



(6) 



Further denote the circulant matrix = circ{0, ■ ■ ■ , 0, 1, 0, ■ ■ ■ , 1, 0, ■ ■ ■ ,0), with corre- 

k times k—1 times 

(k) 

sponding eigenvalues /3j , < i < n — 1. The eigenvalues and eigenvectors for general circulant 
matrices have been well studied, see for example Hoory, Linial, and Wigderson's workP^. Direct 
application yields the eigenvalues and eigenvectors for C^. 



When n is odd, 



A 



(fc) 



2 cos( 



.0 < i < n -1. 



n 



(7) 



When n is even. 



(k) 



2 cos 



' 2nik 



),0 <i <n-l,k 



A 



(k) 



-1)\0 <i<n-l,k 



(8) 



We need to pay attention to the parity of n. By definition, is a circulant matrix generated 
by the vector (0, ■ ■ ■ , 0, 1, 0, ■ ■ ■ , 1, 0, ■ ■ ■ ,0). So the vector that generates has all its entries 



k times 



fc— 1 times 
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equal to zero except the (/c + l)th entry and {n — k + l)th entry, with the value of both entries equal 
to 1. The special case is when n is even and k = n/2, these two entries overlap with each other, 
and the vector has all its entries equal to zero except the (n/2 + l)th entry. 

Now we are ready to solve for the eigenvalue spectrum of R{M, N). We begin with the 
simple example of square rings, i.e., the M = cases. The 2N x 2N adjacency matrix of i?(0, A^) 
has the block form 



A{R{0,N)) 



Cjf In 



(9) 



where In is \he. N x N identity matrix. 



We look for an eigenvector for i?(0, N) of the form Wi = {aivl, vlY where is a constant 
that will be determined later and Vi is an eigenvector of Cj^ corresponding to To be more 
specific, Vi = (1, C\ C^*) ■ ■ ■ ) C^"~^^')S < i < n — 1, where ( is the nth root of unity. If two 
distinct values for a, are found, for any < i < N — 1, then we are done with our search for 
eigenvalues and eigenvectors. 

With such construction for Wi, we need A(depending on i), such that 



Clf In 



\ 



ttiVi 



A 



(10) 
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and so we get the system 



Since ^ 0, otherwise ^^ = , which is a contradiction, so we have 



ai/3f Vj + Vi = XttiVi 



Hence, (aiP^^^ + 1 — a? — Pl'''ai)vi — 0. But then vl Ohy definition, so 1 — a 
conclude that 

a,; = ±1. 



Plug the a\s back into (1 l),we get: 



C'iv 
In Cn 



( m ^ 



+ 1) 



Similarly, 



Cn In 



In Cjf 



\ 



( y 

-Vi 



- - 1) 



/ \ 

— w 
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With suitable labelling, we get that: 

X2i = /3f ^ + 1, \X2i) 



,i = 0,--- ,N -1 



V. 



,i = 0, 



Combining (7) and (8),the eigenvalues of R(0,N) are given by 



2cos(— ) ± 1,0 < i < A^- 1. 



(16) 



(17) 



Now we solve for the eigenvalue spectrum of a general twisted ring R{M, N)(M > 0). The 
2N X 2N adjacency matrix of R{M, N) has the block form 



A{R{M, N)) 



(18) 



We first consider the case ofd — gcd{N, M) — 1. Since d—1, is the adjacency matrix 
of a cycle graph isomorphic to C\,, and so it is similar to C]^. That is there exists a permutation 
matrix P, such that P~^C^ P — Cj^. This implies that the two matrices will have the same eigen- 
values and eigenvectors. Then j3i ' , fi] are eigenvalues corresponding to the same eigenvector, 
say ifi = (1, Cn-> " ' ■> Cw^"^''*)*- We are looking for an eigenvector for A{R{M, N)) of the form 
Wi — {aiVi, Vi)*, where a, will be determined later. If two distinct values for are found, for any 
< i < -/V — 1, we are done with our search for eigenvalues and eigenvectors. 
With this construction for wi, we need A(depending on i), such that 
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a.;V 



( y 

ttiVi 



and so we get the system 



Since vl ^ 0, we have 



But then /3f ^ = can be shown to be included in the cases when 1 — af — 0. So Ui 



Plug the a[s back into (19), we get: 



\ 
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Similarly, 



( \ 



(25) 



With suitable labelling, we get that: 



A2i+1 — — |A2i+l) 



/ \ 

V, 



, i = 0, • • • , TV - 1 



(26) 



The case of d > 1 is treated similarly. The eigenvectors Wi must have the form Wi — 
{aivl', a^vl, ■■■ , Qdvl', vl) with vl as before and vl' = (1, Cn,--- ,Cn^ ~^^*)*> = f > for some 
appropriate multipliers a,. A similar system to do for the d = 1 case will be obtained, and the same 
polynomial whose roots are the eigenvalues Xi will be found. 



Combining (7) and (8), the eigenvalues of R{M, N){M > 0) are given by 



When N is odd. 



2 cos( — ) ± 2 cos(— < i < iV - 1. 



(27) 
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When is even, 



I 



2cos(^) ±2cos( 



2cos(^) ± (-1)\0 < i < N 



2niM 
N 



),0<i<N 



1,M ^ N/2 



1,M = N/2 



(28) 



Construction of Mobius Strips In this appendix, we illustrate how to construct Mobius Strips 
from the respective twisted rings. We use R(2,7) as an example. 



We first remove 6=2+2 x2(and for general R(M,N), 2M+2 ) links. For this particular la- 
belling as in Fig. 6, we first remove the two horizontal links 01 — 07 and 08 — 14, so that the ring 
becomes a band. We then remove 01 - 13 , 02 - 14, 06 - 08, 07 - 09 (and for general R(M,N), 
these are a total of 2M links), so that we can now do the twisting. We then twist one end of the 
band, say the right-hand side. Finally, we rewire back the removed links. And in this particular 
case, 01 - 14, 08 - 07, 01 - 06, 02 - 07, 08 - 13, 09 - 14. The construction for general Mb(M,N) 
is analogues to the above procedure with possibly more links to be removed and glued back. 

Stability and reliability studies of Mobius Strips We give an example to show how the conclu- 
sion above may be applied to general cases. Let's do some rewirings of the network to form Mobius 
strips, which can be regarded as modification of the twisted rings. The number of nodes and links 
and twists are unaffected. We start with a 7-unit rings. The methods for constructing the Mobius 




Figure 6: R(2,7) 
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Strips are shown in supplementary information (II). The stability and reliability landscapes for the 
twisted rings and the respective Mobius strips are shown in Fig.7 and Fig. 8. The spectral gaps are 
shown in Table 1 . Table 2 is a combined ranking of spectral gaps and Stability and Reliability of 
all the 7-unit rings and their respective Mobius Strips. As can be seen from Table 2 as well as 
Fig.7 and Fig. 8, the strict anti-correlation between stability and spectral gap and the strict positive 
correlation between reliability and spectral gap are obeyed not only within the twisted rings and 
Mobius Strips, but also in the combined system. This is a significant extension of the original re- 
sult, as rewirings break the original symmetry of the twisted-ring ensembles: the original C2X 
symmetry is no longer valid in the generalized ensemble, yet our conclusion still holds. Thus, 
we may expect that the conclusion holds for general networks with similar connectance, since the 
heterogeneity in the node and link distribution can always be smoothened by rewiring. A more 
challenging example is the rewiring of a given large network to form a Erdos-Renyi network. We 
find that the usage of the gap ordering for stability and reliability stiU works. 
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Figure 7: Stability landscape of 7-unit rings and Figure 8: Reliability landscape of 7-unit rings 
respective Mobius Strips with a = 0.1. and respective Mobius Strips with p = 0.5. 

17 



Table 1: Spectral Gap of 7-Unit Rings and respective Mobius Strips 
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Table 2: Spectral Gap and Stability & Reliability Rankings 
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